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出 版 前 言 


2011 年 12 月 11 日 是 西安 交通 大 学 杰出 校友 钱学森 先生 
gem. 为 缅怀 钱学森 学 长 ,学 习 他 的 科学 思想 和 卓越 风 
范 , 展 示 其 丰功伟绩 和 人 格 魅力 , 西 安 交通 大 学 举办 了 * 纪 念 钱 
علد‎ ARE 100 周年 "系列 活动 ;作为 制 片 方 之 一 ,参与 西部 电影 
集团 摄制 传记 故事 片 (钱学森 ); 与 中 央 电 视 台 合作 ,出 品 纪录 上 
《实验 班 的 故事 一 沿 着 钱学森 走 过 的 路 ); 扩 建 钱学森 生平 业 
绩 展 馆 ,向 校内 外 开放 ;举办 钱学森 科学 与 教育 思想 研讨 会 ;出 
版 发 行 (钱学森 力学 手稿 )《 钱 学 森 年 谱 ( 初 编 ))《 钱 学 森 第 六 
次 产业 革命 思想 探 微 丛 书 ) 等 。 

钱学森 先生 在 美国 深造 和 工作 期 间 留 下 大 量 珍贵 手稿 ,这 
些 手稿 真实 展示 了 钱学森 先生 博大 精深 的 学 识 .开拓 求实 的 精 
i fu xe x fr r fn fr xi. FE En Er MJ E f r in i f fx in ir ir fh E , 
体现 。 这 里 ,我 们 将 部 分 原稿 整理 汇集 成 册 , 出 版 (钱学森 力学 
FERD AE W RET c cc ANAL. 

《钱学森 力学 手稿 ) 共 10 卷 ,包含 两 部 分 内 容 。 第 一 部 分 是 
草稿 ,包括 扁 壳 , 球 过 和 圆柱 壳 届 曲 分 析 的 公式 推导 和 数值 演 
算 。 在 研究 圆柱 壳 轴 压 届 曲 问题 时 ,为 了 求 得 圆柱 壳 体 的 临界 
压力 ,在 有 关 的 五 百 多 页 章 稿 中 ,对 多 达 二 十 多 种 可 能 的 届 曲 模 


态 逐 一 进行 公式 推演 和 数值 计算 ,最 终 才 找到 满意 的 并 在 论文 
中 采用 的 屈曲 模 态 。 仔 细 观 察 草稿 中 的 数据 列表 ,每 个 数字 有 
效 位 数 都 长 达 八 位 ,在 手 摇 机 械 式 计算 机 作为 主要 计算 工具 的 
年 代 , 这 串 串 数字 凝聚 着 多 少 现今 难以 想象 的 艰辛 劳动 。 

第 二 部 分 是 手稿 ,以 航空 航天 工程 为 核心 ,涵盖 空气 动力 
学 .固体 力学 ,火箭 技术 ,工程 控制 论 和 物理 力学 等 领域 的 部 分 
学 术 论文 手稿 .打印 稿 和 讲义 。 

《钱学森 力学 手稿 ) 是 在 西安 交通 大 学 校 领导 的 大 力 支 持 
FD KAFE U RAE 3 Bt U WHE HUE KE ESE Rk. HB 
出 版 过 程 中 得 到 了 西安 交通 大 学 党 委 宣传 部 、 校 友 关 系 发 展 部 、 
图 书馆 ,航天 航空 学 院 等 的 积极 协助 ,在 此 深 表 感谢 。 
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INTRODUCTION 


I. 


Emissivity calculations for diatomic gases from spectroscopic data 
were developed recently by S. 5. Penner (Ref. 1). His method is based 
upon the use of an average absorption coefficient for the entire funda- 
mental and higher vibration-rotation bands, The method is thus effective 
í when there are extensive overlapping and broadening of the spectral 

lines, and hence is accurate for gases at high total pressures and 

. temperatures, At low pressures, the lines do not overlap and a different 

approach to the problem should be made, Penner and M, H, Ostrander 
I (Ref. 2) have computed the emissivity of carbon monoxide for the case 
of non-overlapping lines by a mmerical procedure, using spectroscopic 
d data obtained recently hy Penner and D. Weber (Ref, 3). The results are 
An excellent agreement with the emissivity determined experimentally by 
W, Ullrich and E, C, Hottel (Ref, 4). The amount of numerical work 
involved is, however, rather heavy. It is the purpose of the present 
paper to develope an anproximate but convenient formula for oaloulafing 
the emissivity of diatomic gases for the case of non-overlapping lines, 


II. FORMULATION OF THE PROBLEM 


If T is the temperature, @ the characteristic temperature, Bee? 
ln 
# the سو‎ mmber, 1 
erptton coefficient at y , p the partial pressure of the emitting gas, | 
| ری‎ wear mande 
y مس‎ 


2 ht 


and L the optical path length, then the emissivity £ of the gas under 
the specified conditions is 
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absorptions for the lines centering on the wave numbers corresponding 
to the indicated transitions, The Se! oan be computed in turn by using 
۸ the results of J, R. Oppenheimer (Ref, 5). As 
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volume per unit pressure, € the effective charge, / the reduced mass, KE EA 
e the velocity of light, Q the complete internal partition 
B's are the internal energy levels given by 


j= 44 l v- * rigoj ge Hg ed] 15) 
En 人- E(0,0) 12 EN = dive Ba I NA Eg moe ETT 


$ 2 De ۵۰ 
ER 4 یڅه‎ rd 


ko 


where x, y » are molecular constante in their standard notations. 
These constants are non-dimensional and are small, The F's and G's are 
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If the fundamental vibration-rotation band gives the main con- 
tribution to the emissivity of the gas, the above equations give the 
necessary information to calculate approximately the emissivity £ 。 


III. APPROXIMATE SOLUTION 


The mmerical work is carrying out the computation indicated in the 


previous section is very heavy. A short formula, however, can be 
developed: First of all, when the lines are به‎ pd from each other, 
sach line can be considered alone, independent of others. Furthermore, 
the value of the factor outside af the bracket in the mmaratar of Eq. (M) 
can be approximated by its value at the center of each line, Thus 
according to 5. S. Penner (Ref, 6) 
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with transitions as indicated, The integrals can be easily evaluated 
(Ref. 6) and are given by the modified Bessel functions I, and Ij: 
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A further approximation can now be made, The magnitude of 
Lee and (۱۶ aro generally quite large if the product pL of pressure 
and optical path length is af the order of unity. Therefore, the 
asymptotic values of the Bessel functions can be used, Then 
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By substituting Eqs. (14) and (15) into (9), the emissivity is calculated 
as a sum over j. 


To carry out the sum over j, one can use the Buler-Maclaurin 
summation formula (Ref, 7), which evaluates the sun by an integral, 
First, due to the smallness of y , $ ,و‎ the following expansions, 
including terms up to the square of D , aro approximate, 
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The corresponding quantities for the transitions j - 1 — j can be very 
easily obtained from Eqs. (16) to ( ty ea 
er this property of symetry, the man of tems from the transition 
3 — 1 - 1 ant the transition j - 1 > 3 for every J is a function of 


T only, Thus, Hager. eee u theme, 
, Ae, Sat Je EP 
- dp ui PM کې‎ E "s Han] 


where 


A is thus a constant independent of temperature and pressure, me fè 
function is simply deduced from the function Q as given by 
Ea. (8): 
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Therefore f is a quantity close to unity, The function g 1s computed 
„FF == 
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The Bian Ad forma مه‎ be now employed to eyeluate 
the mm in the entesivity £ . The resulting integral over j extends 
from 1 to ده‎ , But this range can be made to be frm O to — by 
sinply dedueting the approximate value of the integral from 0 to 1 fran 
the extended integral, Thus 
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The rf) has the mmerical value of 1,225. Finally then the expression 
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where f and g are functions given previously in Eqs, (23) and (24). 
Since the value af f is nearly unity and tho factor before 
g in Eq 2 is small, a good approximate equation for t he emissivity 


£a eho ray HET 
IV. APPLICATION TO CARBON MONOXIDE 


For carbon monoxide, the molecular eonstants are 
Ou ومد‎ 
y= 242.3 کی‎ 
00 Ggs 
4 = 0,0091 
X = 0 
The value of À computed from the neasurenents of Penner and Weber 
(Ref. 3) 1s 22 
A * E مې‎ 
They have also determined (Ref. d) > to be 0.077 an" at ono mtnosphero 
of total pressure, 


According to the aprroxinate equation a, the emissivity at 
T = 300% and a total pressure of ons atmosphere 


۱.۸۱0۵۳ 


0 
1.092 + 
£ = we “مدع‎ det (29) 
where pL is in atm-cm, By using the nore exact equation (25), the 


emissivity is 
Bus 2 
ع‎ m 60 
Mo 
The difference between the anproximate value and the more exact value is 


quite small, The comparison between the comouted emissivity and the 
measurements of Ullrich and Hottel (Ref. 4) is show, in Fig. 1, The 
agreement is quite satisfactory up to pL of approximately 10, 
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This paper shovs that the combustion in the rocket motor can be stabilized , 
against any value of time lag in combustion by a feedback servo link from a 
chamber pressure pickup, through an appropriately designed amplifier amd & 4 « 
control capacitance, meting on the propellant feediine. The technique of 
stability analysis 1s based uyon a combination که‎ fiche diagram and an 
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The phenomenon of rough burning in liquid propellant rocket motor has 


been interpreted as the instability of the coupled aysten of propellant feed 
and combustion chamber by D. F. Ounder and D. R. Friant (Bef. 1), H. Yachter 
(Ref. 2), M. Summerfield (Ref. 3), and L. Crocco (Ref. 4). The essential 
feature of these theories is the time lag between the instant of injection of 
the propellant and the Instant when the propellant is bumed into hot gas. 
Orocco has further improved on this concept by considering the tine lag as an 


integrated effect of consecutive stages, each of which is controlled by the 
prevailing pressure in the combustion chamber. As a result of this new concept, 
Orocco showed the possibility of intrinsic instability with constant injection 
rate not influenced ty the chamber pressure. 

The present paper will first give a slightly more general formulation of 
Oroeco's concept of time lag, allowing arbitrary pressure dependence of lag. 
Then the problem of intrinsic stability is Aiscussed by applying a mothod 
sugrestod by M. Satche (Ref. 5). This method is based upon a modification of 
the Nyquist diagram and is particularly useful for systems having time lag. Tor 
easy reference, this now diagram vill be called the Satche diagram. The later 
sections of the paper will show the possibility of stabilizing the combustion 
by moans of a faodhaor servo for all values of tine Ing. Such possibility of 
servo-stebillsation was first mentioned by Y. Bollay in his admirable paper 
(Raf. 6) on the application of servo-mechanisus to aeronautics. The present 
study definitely shows the power of this idoa. 


Time Laz in Conbustion 

Let m(t) ba the mass rate of generation of hot gas by combustion at tine 
instant $. Consider, for simplicity, a monopropollant motor. ‘Thon the mass rate 
of injection at t can be denoted by tt). Let (8) be the timo lag for that 


pareel of propellant which is burned at the instant t. Then the mass burned 
during the interval from $ to $ + dt must be equal to the mass injected during 
the time from t= tot- +d(t= ). Thus 


a) 


The mass of hot gas generated is either used to fill the combustion chamber 
dy raising its pressure p(t) or is discharged through the rockot nozzle. If the 
frequency of the possible oscillations in the chamber ie small, then the pressure 
in the chamber can be considered as uniform and as a first approximation (Ref. 7) 
the rate of flow through the nozzle can be taken as proportional to the 
instantaneous chamber pressure p(t). Thus if is tho stendy mass rate flow 
through the system, is the average mass of hot gas in the chamber, and if 
the volume occupied by the unburned liquid propellant is neglected, 


(2) 


where ie the steady state pressure in the combustion chamber, 
By following Crocco, the non-dimensional variables for the chamber pressure 
and the rate of injection are defined as 


(3) 


and are then the fractional deviation of pressure and injection rate 
from the average. With Eq. (3), can be eliminated from Eqs. (1) and (2), 
and 


(4) 


fo calculate the quantity „ 0740٥0 concept of pressure dependence 
of tino lag has to be introduced. If the rate at which the liquid propellant is 


— 


prepared for the final rapid transformation into hot gas is a funotion 
then the lag 1s determined by 


= Constant 


By differentiating Eq. (5) with respect to t, 


The concept of smell porturbation from the steady state will now be explicitly 
introduced: assumo that tho deviation of the pressure p fron the steady state 
value is small. Then at the instant t and at the instant 

te can be expanded as Taylor's series around + By taking only the first 
order terus, 


Here 48 the lag at the average pressure + a constant now. Then 


By combining Zqs. (4) and (6), the folloving equation ia obtained 


(9) 
If n is a constant independent of + then is proportional to ۰ 
This is the fom of assumed by Crocco. The present formulation of the 
problem is slightly more general in that . is arbitrary and the value of n 
1s to be computed by using Eq. (8), and is a function of ۰ is, of course, 
the gas transit time. 


Croooe called the instability of combustion with constant rate of injection 
the intrinsic instability, If the injection rate is constant not influenced ty 
the chamber pressure p, then ۰ EF A 
dy the following 的 nple equation obtained from Nq. (7), 


Thea an 


This is the equation for the exponent a. 

Crocco determined the value of the complex number a by studying the set of 
two equations for the real and the imaginary parte of Eq. (11). However if the 
point of interest is whether the system is stable or not, one can use the well- 
known Cauchy's theorem vith advantage. Let 


6 


Then the question of stability is determined by whether has zeros in 
the right half of the complex e-plane. This question itself can be in turn 
answered by watching the argument of when s trace a contour enclosing 
the right half s-plene. Specifically, let s trace clockwisely the contour 
consisted of the imaginary axis and a large half circle to the right of the 
imeginary axis (Fig. 1). If the vector make a number of complete 
clockwise revolutions, then that number is, according to Cauchy's thoorem, the 
difference between the number of zeros and the number of poles of in the 
Tight half seplane, Since evidently has no poles in the seplane, the 
number of revolutions of is the number of zeros, Hence for stability, 
the vector must not make any complete revolutions, as s traces the 
specified contour. Therefore the stability question can be answered by plotting 
graphically on the complex plane. This graph is, of course, the well 
known Nyquist disgrem. 

A direct application of this method to given by Eq. (12) is however 
inconvenient for th: @omplication caused hy lag term (Raf. 8). Me Satche 
(Ref. 5) however proposed a vory elegant and ingenious method to treat such 


system with tine lagi Instead of + break 1 into two parts, 

(13) 
vhere 

Q4) 
The veotor is thus a vector with vortex in and ito tail on ۰ 


The graph of is the unit circle for s on the imaginary exis. For s on the 
large half circle, 15 within the unit circle, The graph of is 


tho straight line (Pig. 2) paralled to the imaginary axis when s is on the 
imaginary axis. When s is on the large half circle, is a half of 
a great circle closing the contour on the left. A momente! reflection will 
ahow that in order for the vector not to make complete revolutions 

for any value of „ the contour must lie completely out of the 

contour, That is, for unconditional intrinsic stability, 


or (15) 


men „ the contour and the contour intersoat. Stability 


is still possible however, if for within the unit circle (Mg. 3), is 
to the right of + This condition is satisfied if 


(17) 
+ ‘Therefore when + has the oscillatory solution with the 
angular frequency *, 

These results on intrinsic stability were obtained by Crocco. The present 
discussion with the Satche diagram however seems to be simpler. Tor the more 
complicated stability problem treated below wm feed-aystem and servo-control, 
the solution is hardly prectieal without the Satche diagram. 


Araton Dymanica vith Serve-oontrak 
Consider now a system including the propellant feed and a servo=control 
represented by Mg. 4. In order to approximate the elasticity of the feed line, 


m 


a spring load capacitance is put at the midway point between the propellant 
pump and the injector. The spring constant is to be computed from the feed 
line dimensions.* Near the injector there is another capacitance controlled 

by the servo. The servo receives its signal from the chamber pressure 

pickup through en amplifier. If the feed system and the motor design is fixed 
by the designer, the question 1s whether it is possible to design an appropriate 
amplifier so that the whole system will be stable. Because there is no accurate 
information on the time lag of combustion, a practical design should specify 


unconditional stability, 1.e., stability for any value of . 


Let be the instantaneous mass flow rate out of the propellant pump 
and be the instantaneous pressure at the outlet of pump. The average flow 
rate must be + The average pressure is + The pomp characteristics can 
be reprosented by the following equation, 


(18) 


If the time rate of change of mass flow is small, ie simply related to the 
slope of the head-volume curve of the pump at constant speed near the standye 
state operating point. For constant pressure pump or the simple pressure feed, 

is zero. For conventional centrifugal pumps, is approximately 1. For 
displacement pumps, is very large. 

Let be the instantaneous mass rate of flow after the spring loaded 
capacitance, be the spring constant of the capacitance and the 
instantaneous pressure at the capacitance. Them 


Qs) 


* See the Appendix for details. 


where is the density of the propellant, a constant. 

In the following calculation, the pressure drop in the line by frictional 
forces vill ve neglected. Then the pressure difference ia due to the 
acceleration of the flow only. That is 


(20) 
where A is the cross-sectional area of the feed line, a constant, and is the 
total length of the feed line. Sinilarly, if de the instantaneous pressure 
at the control capacitance, 

(a) 
If the mass capacity of the control capacitance is C, then 

(22) 
Since the control capacitance is very close to the injoctor, the inertia of the 


mass of propellant between the control capacitance and the injector is 
negligitle. Then 


(es) 


where Ay tho effective orifice area of the injector. يد‎ can be elininated 
from the calculation hy noting that at steady state, tho difforonoe of pressures 
and „or is 


(24) 


Zquations (18) to (24) describe the dynamice of the food system, Ay a 
straightforward process of olimination of variables, a relation between 


end 7 is obtained. To express this relation in non-dimensional form, the 


folloving quantities are introduced, following the notation of Oroceo: 


(26) 


and (26) 


where the gas transit tine given dy Ba. (9). Thon the non-dimensional 
equation relating and HI 


here is the non-dimensional time variable defined ty Ba. (9). 

The dynamics of the servo-control is specified by the composite of the 
instrument characteristics of the pressure pickup, the response of the amplifier 
and the properties of the servo. Since it is not the purpose of the present 
paper to discuss the detailed design of the servo-control, the overall dynamics 
of the sorvo-control is represented by the following operator equations 


(28) 


whore Y is tho ratio of two polynomials with the denominator of higher order 
than the numorator. 

Equations (7), (27) and (28) are the three equations for the three 
variables ۰ and . Since they are equations with constant coefficients, 
the appropriate forms for the variables are 


(23) 


By substituting zq. (29) into Equ. (7), (27) and (28), three homogeneous 
equations for a, b and c are obtained. In order for a, b, c to be non-zero, 
the determinant fomed by their coefficients must vanish. This condition can 
be written ag follows: 


This is the oquation for determining the exponent s. F(s) is now recognized as 
the overall transfer function of the servo-ooritrel link. The complete system 


stability depends upon whether Eq. (3) gives roots that have positive real part. 


Instability vithous Servo-Controk 

The systom characteristics without the servo-control can be simply obtained 
from the basie equation (30) ty setting F(6) = 0. Let it be assumed that the 
polynomial multiplied inte has no zero in the positive half s-plane, as is 
usually the case. Then 50. (30) can be divided by that polynominal vithout 
introducing poles in the positive half e-plane into the resultant function. 
That is, for the Satche diagran, one has again 


is thus egain the “unit circle". 1s nov mach more complicated: 


a) 


when s is pare imaginary is given by setting s= O in 


The intercept of 
مولا‎ (831), 4 


(aa) 


Since all the parameters ۰ » P positive, the magnitude of is 
now smaller than the magnitude of given by Zq. (14) for the intrinsic 
stability problem, Thus the effect of the food-aysten is to move the 
curro tovards the unit circle of in the Satche diagram. For instance, for 

۰ is just tangent to the unit circle for the intrinsic system 
without considering the propellant feed. But vith the propellant food-ayatem, 

contour will intersect the unit circle and the system will become unstable 
for tino lag — exceed a certain finite value. The influence of the feod-ayaten 4 
is thus always de-stabilising. Thie is further confirmed by considering the 4 
asymptote of for large imaginary — . obtained from Ba. (31), That is 


(33) ۳ 


Therefore for large inaginary a, approaches asymptotically a line 
parallel to the imaginary arie at a distance 


to the left of the imaginary axis. The effect of fesd-aystem is again to nove 
towards the unit circle. 


It is thus evident that for the parameter mear 1/2 or larger than 1/2, 
4% would be impossible to design the system for unconditional stability, In 
the Satche disgran, contour and will always intersect without a 
servo-control. 


Somplete Stability with Servo-Gentro) 


If the polynomial E(s), 


" which multiplies into in Eq. (30), has no poles and seros in the right 
half e-plane, then the ocouranse seros of the expression in 24. (30) in the 
Tight half s-plans can be determined from the Satche diagram vith 


(35) 


As a traces the contour of Pig. l, is again a unit circle. Therefore 
| 17 simultaneously the contour is completely outside the unit circle, 

١ there can be no root of Eq. (30) in the right half s-plane, In other words, if 
| the transfer function P(s) of the servo-control link is so designed as to place 
1 the contour completely out of the anit circle (Fig. 5), then the system 
۱ is stabilised for all tine lags. 


Of prinary interest is the behavior of when s is a pure imaginary number 
۰ real. Thus 


‘This contour for is plotted in Fig. 6. It is evident that for sufficiently 
large values of tino lag, the system will be unstable. On the other hand, if 
the contour can be changed ty the servo-control to say 


Then as plotted in Fig. 6, the new contour is complotely outeide of the 
unit oirole of + Therefore the system is now unconditionally stable. A 
straightforvard calculation from Bag. (31) and (35), shows that the required 
transfer function F(s) for the serve link is 


The servo link has thus the character of an integrating circuit. If with given 
response of the chamber pressure piekup and of the servo for the control 


capacitance, an amplifier could be designed to give an overall transfer function 
close to that specified above, the conbustion can de stabilised by such a servo- 
control. 

As the second example, taks 


Since s the feed pressure is thus constant with evon variable flow 
of propellant. The case then corresponds to that of a simple pressure feed, 
Without the servo-control, 


This contour of 1s plotted in Pig. 7. It is evident that without servo= 
control the eonbustion will be unstable for sufficiently long time lag. In fact, 
the stan is even loss stable than the system considered in the first ardet 
It will become unstable at shorter time lag. The part of the contour near 


2 is of special interest. Near „ the contour comes so close of 


the unit circle of that 1۶ the value of tine lag is such aa to make 
and for 2 very close to each other, then an almost undanped oscillation 
at 2 can occur. This critical valus of is evidently smaller than 
the critionl determined from the true intersection of with the unit 


círole at 0.65. Such near instability at smaller values of tino lag 
هده‎ be easily ovarlooked in the analytic treatment of the stability condition 


by Crocco, and yet such possible instability should not be dismissed. This, 
perhaps, indicates the supariority of the present graphical method. 
should be displaced out of the unit 


The roqui red 


For unconditional stability, 


circle, to, say, the sano "stable" contour as in the first exemple, 


transfer function F(s) {o calculated to be 


The required servo link must then have the character of double integrating 
circuit. Furthermore, the transfer function has two purely imaginary poles at 

+ This unreslistic requirement on the amplifier conos from the original 
feod-systen dynamics and is due to the neglection of frictional damping in the 
feed line. In any actual system, the frictional damping in the feed line will 
remove these purely inaginary poles of the required transfer function F(s) and 
replace them by two complex conjugate poles. 


In the preceding discussion of servo=stabilization, it is assumed that the 
polynomial H(s), Eq. (34) has no pole or sero in the richt half seplane. This 

is however not necessarily the case. In general then, one should first investigate 
the number of zeros and poles of H(s) in the right half seplane. To do this, it 
should be recognized that the polynomial in 24. (34) before the factor (و)7‎ 
usually does not have dico dx de right half s-plane. Therefore instead of 
studying H(e), one can study the ratio of H(s) and that polynomial. That is, 


the number of zeros and poles of E(s) in the right half s-plane is the sane 
es the number of zeros and poles of the following function 


(37) 


According to the Nyquist criterion, the number of poles and zeros for leX(s) 
in the right half plane can be found by plotting the Nyquist diagram of leX(s) 
vith s tracing the contour of Fig. 1. In fact, if 14K(s) or Als) has zeros 
and poles in right half s-plano, than K(s) will carry out clockwise 
revolutions around the point «l, ns s traces tho contour of Fig. l. Range the 
necessary information on Bis) ean be obtained hy plotting the Nyquist diagram of 
0 x(a). : 
When ono dividos the Eq. (30) by Rte) in order to obtain and as 
given by Eq. (35), zeros and poles are introduced in the right half 
‘plane. The poles of K(s) mst come from F(s), since the polynomial in the 
denominator of Zq. (37) has no sere in the right half seplane. Therefore the 
original expression in 24. (30) also has poles in the right half s-plane, 
Hence in order for the original expression in 2q. (30) to have no zero in tha 
right half plane, must make clockwise revolutions 
E around the unit circle. In order the stability be unconditional, 1.0,, stable 


for all tine lag, the contour should never intersect the unit circle. 
Therefore the general unconditional stability criteria are, first, 


contour 
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completely outside of the unit circle; and, second, making counter- 
clockwise revolutions around the unit circle as s traces the conventional 
contour enclosing the right half s-plane, These are the criteria for stability 
with the Satche diagram. To determine , one has to use the Nyquist diagram 
of K(s), Eq. (37). Thus the stability problem for the genoral case roguiree 


both the Satche diagrea and the Nyquist diagrams. (Fig. 8) 


Goncludine Renarka 

In the previous sections of this paper, the theoretical possibility of 
completely stabilizing the combustion for any value of time lag by servo=control 
is demonstrated. The great flexibility of electronic amplifier seons to 
indicate that this thooretical possibility cen be always realized. On the other 
hand, without the servo link, unconditional stability is shown to be generally 
impossible. Therefore the concept of feedback servo is indeed a powerful tool 
in controlling the behavior of a time-lag system, It 13 to be realized, of 
course, that the proposed scheme is but one among many. No attempt is made here 
to give an exhaust treatment of all possibla schemes. The best scheme is 
certainly to be determined by detailed considerations on all aspects of the 
engineering problem, such as the possibility of high frequency accoustic 
oscillations which are not considered here. The main purpose here is to give a 
general discussion of the concept together with a suggested general method of 
analysing the stability by the Satche diagram. 

It is of interest to point out that stabilization by servo-control is only 
one phase of the general concept of feedback link. The opposite case of 
de=atabiligation could be of importance also, For instance, consider the 
so-called velveless pulsejet. It is not always possible tc operate the engine 
with the desired pulsation. With a feedback servo linking the combustion chamber 


we 
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pressure pickup through an amplifier to the fuel line, the system can be de- 

stebilized at the dest red operating frequency and thus operste the ongine at 

that frequency of pulsation. his application of servo=destabilization gives 
the vulveless pulsejet a new flexibility and an extended range of operation. 

Therefore 1% moens worthwhile to explore carefully all possible applications 

of feedback control to systems with tine lag. 


Salenlation of Parametere تر‎ end 3 


If Le and e" are the characteristic length and the characteristic velocity 
of the motor, and if Ta is the chauber temperature, R the gas constant, the 
transit time ie 


To calculate J and 3 defined by Bq. (25), it is more convenient to use the 
averege propellant velocity in the feed line. hus 


Thus, according to Eq. (25) 


A consistant set of units would be An slugs per cubic foot, in feet 
por second, in feet, in seoonds and in pounds por square foot. 
If d tho diameter of the feed line, h ite thickness and 25 the Youngs 
modulus of the tube material, then „ the change in volume of the food line 

per unit rise in pressure, is 


Therefore By. (25) gives 


A consistant set of units would de ín pounds per square inchs Z! in 
pounds per square inch, 1n feet, in seconds, and in feet per second. 
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Figure Captions 


Contour traced by the variable s for the Satohe Disgran or the 
Nyquist Diagram. 


Stable Satche Diagram for Intrinsic Oscillations; 

Unstable Setche Diagram for Intrinsic Oscillations; 
Servo-controlled Liquid Monepropellant Rocket Motor. 

Satche Diagram for the Original end for the Servo+Stadilized System, 
Satche 7 lee sa وي‎ OE وو‎ 


P=3/2, Jwa B= 1/4 -1 


without servo intersecte the unit circle; with servo 
io outside the unit cirele, unbors beside points are the value of 


Satche Diagram for the Original and for the Sorvo-Stabilized System 
P=3/2, J= 4, B= 1/4, = 0 


without servo intersects the unit circle; with servo 
is outside the unit circle. Numbers beside points are the value of 


curve for jive dotted curve for negative 
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Section 3 


Engineering Cybernetics 


Conflicionts.‏ تس 
Let us consider the simplest system - a firat order system,‏ 
That is, the differential equation of the system is a firat order linear‏ 
differential equation of constant coefficients, If the system is assumed‏ 
to be free and is not subjected to "forcing function then the differential‏ 
equation can be written as‏ 


„ Ou 


4 may be called the spring constant and is real. When there is no 
variation of & with respect to time, It vanishes and then Eq. (1.1) 
requires 4=0. Therefore the stationary state or the equilibrium state 
of the system corresponds to gwo. 

The solution of Eg. (1.1) ie 


و 


W R or 
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Y” Y io thus the initial disturbance of the system from the equilibrium ue 
state. The behavior of the system for (%, is illustrated in Fig. 1.1 
for both positive É and for negative £ . It is seen that for 4% 
the magnitude of Y decreases with time. Then as time increases 
indefinitely, 多 -0 . “Therefore 


for Coe, the disturbance of the system will eventually disappear. 
The system can then be said to be stable. When {<o , the disturbed 
motion of the system increases with time and eventually the disturb- 
ance will become very large no matter how small the initial displace- 
ment is, and will never return to the equilibrium state once disturbed. 
Such systems are thus unstable. 

For systems of higher order, the differential equation will have 
higher derivatives. The n-th order system has the differential equation 


3 + =o as 


For a physical system, the coefficients يچ‎ ٠۰۰ كيه‎ are real. Then 
the solution of Eq. (1.4) can be written as 


^ 
aer a.) 
where o4 , fp are real and are related to the coefficients Go,» de 
and fp are the phase angles. The motion of the system is thus only 
stable if all Ku e negative. If one of them is positive, the disturb- 
ance will eventually diverge, and the system is thus unstable. 

From the above examples it is seen that the crucial question to 
ask about the behavior of a linear system of constant coefficients is the 
question of stability. Needless to say, the usual aim of an engineering 
design is stability. The question of stability can be answered however , 
once the coefficients of the differential equation are specified. Incase 
of the simple first order system specified by Eq. (1.1), the only infor- 
mation that matters is the sign of the coefficient K. 


t> Linear System with Variable Coefficients 
H there is a variable parameter in the system under study, the 
stationary or the equilibrium state of the system can be changed by 
changing this parameter. It is natural then to expect the coefficients of 
the linear differential equation describing system to be also functions 
of this parameter. For instance, the aerodynamic forces acting on an 
aircraft are functions of the speed of the aircraft. If the speed of the 
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Chapter II 
Method of Laplace Transform 


For linear differential equations with constant coefficients and 
with time ځ‎ as the independent variable, the method of Laplace trans- 
form is particularly useful in finding the solution. Of course, the 
problem can be solved by a number of other methods; but the method of 
Laplace transform appeals specially to the engineering scientists in that 
it reduces all problems to a uniform basis. The procedure of solution is 
then standardized and a general approach is possible. The theory and 
practice of Laplace transform is discussed in many texts. * It is not the 
purpose of the present chapter to do this. The purpose here is rather to 
give a summary of results which are useful to our discussion in the sub- 
sequent chapters for easy reference. For details andproofs, the reader 
should consult the texts cited. 


A. Laplace Transform and Inversion Formula 
If lc) is a function of time variable É defined for £>0 then 


the Laplace transform Y(S) of YO) هذ‎ defined anto 
YG) Je? que) dé (2.1) 
6 


where 5 is a complex variable having a positive real part, 2>0 
For otber values of 5 , the function le) is defined by the analytic 
continuation. The dimension of Y($) is the dimension of Y multi- 
plied by time. 

When Y(S) is known, the original function for which Y(s) is the 
Laplace transform can be obtained in all cases by the Inversion Formula: 


* See for instance, H. S. Carslaw and J. C. J. 
Methods in Applied Mathematics", Oxford,(194! . 
"Modern Operational Methods in Engineering". McGraw Hill. (1944). 
For morgcomplete theory, one should consult G. Doetsch, "Theorie 
und Anwe der Laplace-Transformation", J. Springer, Berlin. 


(1937); or D. V. Widder, "The Laplace Transform", Princeton, (1946). 


## We shall use throughout capital alphabet to denote the Laplace trans- 
form of quantities denoted by a lower case alphabet. 


Therefore the error signal vanishes as - 
Consider now another example of the input: Let the input be 


sinusoidal, la 
xlt) = xg, e` 


where 2 is the amplitude and w is the frequency. Then 


XO) 3 (3.12) 


The output due to the initial condition is the same as before. The output 
due to input is given by 


4 
Ort TEE ولېه‎ +1] 
Therefore according to our dictionary, the output yo is 
Y 
+ E E 
MEES EE S 


The first term is a pure subsidence and the second term is the steady 
state output. Thus سو‎ ca 


Del, T Z = Fü) 


This is in full agreement with our general result given in Eq. (2. 16). 


Since 


dag CU (3.13) 
— 9 
the steady state output can be expressed as 
om] 


[ze] > n L 


Therefore the amplitude of the steady state output is reduced by the eg 
factor , + ae) in comparison with the input, and the phase of * 
the output lage behind the input by the amount dan ( (w5). For low 
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[2] 
the aileron deflection $. equation for the roll angle Pm 


E 


Now let [TE Ae then the above equation becomes 


48 = فیط + فك + 


. H the roll speed is zero at ¢=0 , then the transformed equation is 
st د-09‎ 0 


The transfer function cs) is thus 


Res £ £ / 
ag Te y repe pm 


The behavior of the system is determined by the transfer function is 
thus similar to the cantilever spring with dashpot and the simple lag 
network. Here the characteristic time T, is MOU If the damping 
is very small, then 5 —=o and the behavior of the system becomes that 
of the simple integrator. 


3.4, Second Order Systems 

En Let us return to the cantilever spring with a dashpot, Fig. 3.1. 
Only now we attach a mass m to the dashpot end, The mass will 
introduce an inertia force m d , and the equation of motion is now 


m 4 we oky bz 
with the initial conditions 
4 


A; Fi (3.37) 


The differential equation of motion can be rewritten in more convenient 
form by introducing the following parameters: 


with G(s) given. The method of Evans determines such roots as 
functions of the gain < . and is thus called the root-locus method. 
When this is done, any set of specifications on the roots gives a proper 
choice of the magnitude of X . This method then goes much beyond 
the mere satisfaction of criterion a) of Sectionf2, but actually solves 
the design problem for all three criteria stated in that section. 

Now let Gs) be specified by ite zeros ff , p. * and 
ite poles . g- f. Them from tha definition of gain given 
by Eqs. (3.16), (3-21) and (3. 23). 


Gc) A اسو‎ BOETE. 


SA] Fa) Wee A 


For physical systems, the polynomials in the numerator and the 
denominator of Ge) has real coefficiente. Then the pa are eitber 

real or form complex conjugate pairs. Similarly the g's are either 
Ale A A مني‎ ed 


mb A wek‏ „ دص A < u‏ /مسم Ph‏ ۴م 


real or form complex conjugate pairs. 


SR لہ‎ bd abo kaihoe. T tha, we shall غفا‎ A te 

A 4 ad Anti. Gal thee dk f Gis) A f بس موه‎ 

biplane verde thane Het f the be, د‎ SÈ m LA ma. afr تیه‎ 
mt the . $. Ep OR) مخ‎ mede di PY 


1 SEN — 27 > N 


کیچ 5-2 
)4.17( 


tos. Th ¿8 
the g tos. le vector © e 


$ is the variable point in the c plane 
By using Equ. (4.18) and (4.19), G(s) can be ne = : 


Gis) = A LA. 2 és 
(2 e“y(Q,e' 


Hence we can write Eq. 


R = AFR” SA e. Gq) (4.21) 


and © z(grg t په )-(ل و‎ ----18,) 14. 22) 


Since P's and Q's are magnitudes of vectors defined by Eqs. (4. 18) 
and (4.19), they are positive, Therefore 2 is positive. The basic 
equation for the roots of inverse system transfer function, Eq. (4, 15). 


E 


KR 
Therefore to satisfy this equation, we must have 


=-/ 


KR =/ 14.23) 
and Oath (4.24) 
” The method of Evans consists of two steps: The first step is to 


determine all 5 that satisfy the appropriate angle condition of Eq. (4. 24). 
Then knowing such root-locus, we can compute & and hence & , by 
Eq. (4.23), for each point on the root-locus. Evans has developed a 
number of useful rules for plotting the root-locus, We shall explain 
those rules presently. 

Rule 1 For X =0, Eq. (4.15) shows that lse. Thus for K . 
the roots of LE) are poles of G(s) or the root-locus starts at the 
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oscillating servomechanisms are less flexible than are oscillating L 

control servomechnisms in which the oscillation is supplied by an : 

independent generator. ۲ 
An elementary precaution to be observed, in order that the ۳ 

curve, which is constrained to pans through the point - 2 , shall -void 

the neighborhood of the point -/ , is that the curve shou).! interest the 

real axis, at the point -2 , perpendicularly. This irn Ales thatthe vector 

1 veel should be varying slowly in amplitude, and capidly in angle, at 
the frequency at which the system oscillates. 


A relay omaclimites is a non/ linear device, But by mixing the 
signal with a sinusoidal oscillation df high frequency and large anplitude, we 
the output is made to be linear with resp:ct to the signal. Thus the 
essential concept of oscillating control se: echanisms is the lineariz- 
ation of a non-linear system. J. M. * has shown that this concept is 
applicable to any non-linear system, anc he calls this method the general 
linearizing process for non-linear contro. systems. We shall call the 
resulting servomechanism the general oscillating control servomechanism. 
Let us consider a general function g(x) where y is the output 
and = is the input. If instead of the varlıble x , we substitute the 
sum 76 where € is much smaller than عد‎ . Then if the function yon 
is regular, we can expand glare into a Taylor series as 


gere سړ-‎ E + ($) K d (6.15) 


We now specify the input x as a periodic fün-tion of time ځ‎ with the ay 
period 7 , and € as a constant. Thef it is clear that gen is also 
a periodic function of time with the samd perod 7 . Sòme is true for 
Bite ana dK" . Periodid fur-tions can be expanded into 
Fourier neries; thus if we neglect powerènf € higher than first, we have 


eo 
Ars a „(, Gs not + يبأ‎ Sim nwt ) 


Mal مه‎ 
+ © ] په‎ +3! (8, Cas mat + b, Sin awé)] (6. 16) 
Mel 


* J. M. Loeb, Annilles des Télécommanications, 5:65-71 (1950). 


We proceed to evaluate the right-hand member of (7. 11) by the method 
of residues. 

The integrand has certain poles, the poles of 4(s) . tying to the 
left of the path of integration, and other poles, which are the roots of the 
equation /- FÉ 8-8), , tying to the right of the path of integration. 
is easily seen that the integration upward along the line I -¿ to, E 
is equivalent to integration in the clockwise direction ub the closed 
contour formed by that line and the infinite semicircle in the right balf- 
plane. Hence the right-hand member of Eq. (7. 11) is -4 times the 
sum of the residues of the integrand with respect to the several roots of 
the equation /- — o. 

Now the typical root of the equation is gs eu where m 
is an integer, and the residue of the intagrand with respect to that pole 
is geer animo . Therefore finally 

LI 


h => | (sertir) : رو ددم‎ as 


* co 


4 
This formula gives considerable insight into the properties of /£($) ۰ 
and at times may be useful in making approximate calculations. How- 
ever, we can easily obtain an exact representation of g in finite 
form. 


The function g lo) can be represented as the sum of a finite 
number of partial fractions. = 


— = ۷ `... 


Chapter VEE 


Linear Systems with Time Lag 


In this chapter. we shall introduce another new element into our 
linear systems with constant coefficient: the time lag. By time lag, Ç 9 
we mean that the relation between the different variables of the system 
cannot be expressed as a relation of these variables all taken at some 
time instant # ; but on the contrary, the relation involves variables. 

some taken at the time instant ¢ , and some taken at an earlier instant 
r Those taken at the instant ¢-Z then lag by the interval C behind 
the variables taken at the instant # . This time lag is thus quite different 
from the characteristic time constant of a first order linear system 
introduced in Section 3.1. Time lag systems are represented by 
differential-difference equations of constant coefficients and are more 
complex than the linear systems studied previously which are represented 
by differential equations. Systems with time lag were studied by many 
investigators; for instance, A. Callander, D. Hartree, and A. Porter* 
and N. Minorsky. ** Our interest here is, however, somewhat more 
restricted, We wish to know: How can we analyse the performance of è 
feedback servomechanism if there is a characteristic time lag Ç in the 
system? We wish, specifically. to modify the method of Nyquist of 
Section 4. 3, to time lag systems. 

We shall developf tbe theory by treating a particular example of 
such systems; the example of stabilizing the combustion in a rocket motor 
by feedback control. The problem of combustion instability in rocket 
motors was treated by many authors, the following analysis of combustion 
lag time originates from the work of L. Crocco. ### For simplicity of 
calculation, **** we shall consider only the case of دد‎ called low frequency 
oscillation in a rocket motor using single liquid propellant. 


* A. Callander, D. Hartree and A. Porter, Phil. Trans. Royal Society 
of London (A), 235:415-444 (1935). 


* N. Minorsky, J. Appl. Mechanics (ASME) 9:67-71 (1942). 
*** L. Crocco, J. American Rocket Society, 21: 163-178 (1951). 


The following discussion is based upon a paper in J. Americas Rocket 
Society, 22:256-262 (1952). 
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Chapter De 
Linear Systems with Stationary Random Inputs 


In the previous chapters, the inputs to a system are considered 
to be definitely specified functions of time ۶ . However, there 
many engineering problems for linear systems with constant coefficients 
where the inputs cannot be so definitely described. An example of such 
engineering problem is the problem of the motion and the stresses 
induced in the structure of an airplane wing in turbulent air stream. 
Here the input can be considered to be the time varying air-flow pattern. 
But the airflow pattern cannot be described as a definite function of time, 
N but has to be recognized as a random function of time, specified by 
certain statistical characteristics. It is then evident, the output of the 

' system, the stresses in this case, must be also a random function and 

can be described also only in statistical terms. The first objective of 

this chapter is then to find a convenient method of calculating the statis- 

tical properties of the output from the specified statistical properties of 

the input. This forms an easy extention of the early investigations by 

P. Langevin of the Brownian motion. 

Another example of random input is the so-called noise in control ¿A£ 
signals, The noise is introduced by the disturbances and the fluctuations 
beyond the control of the designer. The problem of noise is a problem 
of much research in connection communications engineering. There 
the cÄntröl question is how to System ef-eede so that the effects 


of the unavoidable noise can be minimized and the us: information of 
discus e 0 r te. 
the نا ناه‎ oot یک‎ 


however, somewhat different: In our problem, the random output is the 
only output of the system. Our parpose in the design of the system, 
particularly the design of (be feedback servomechanism. is to obtain with 
A given input an output of the desired statistical characteristics. We shall 
see that the method of transfer function developed in the previous chapters 
remain useful in the present task. 


٧ Statistical Description of د‎ Random Function 


Let us consider a system which generates a random function yo. 
Now to formulate the concept of š statistical description of such a random 


Chapter 1 
Linear System with Variable Coefficients 


The only system with time varying coefficients considered in 
detail in the previous chapters is the pendulum with a periodic force at 
the supporting end, discussed in connection with the phenomenon of 
parametric excitation and damping. All other systems considered do 
not have coefficients of their differential equations that are explicitly 
functions of time. We bave, however, shown in Chapter f, that linear 
systems with time varying coefficients can have behavior entirely different 
from systems with constant coefficients. In this chapter, we shall again 
take up this question and discuss in some detail such a typical but simple 
system: the short range artillery rocket. We shall demonstrate that the 
question of stability of such a system with variable coefficients cannot be 
solved in the same manner as the linear system with constant coefficients. 
Not only is the method of Laplace transform and transfer function useless 
for the purpose, but also we are forced to change our entire approach to 
the problem. 

We shall study the motion of a fin-stabilized artillery rocket 
during the period of action of the rocket thrust. We shall be particularly 
concerned about the angular deviations of the rocket axis from the launch- 
ing angle due to the action of disturbances when leaving the launcher and 
the subsequent damping action of the fins. The general problem of 
dynamics of artillery rockets has been studied in great detail by various 
authors in different countries during World War U. The American work 
is summarized by Rosser, Newton, ani Gros. * The work done in 
England is reported by Rankin, ** Carriere's paper*9* represents 
French investigation on the same subject. Our discussion here 
will be greatly simplified and has the purpose of only bringing out 


۰ J. B. Rosser, R. R. Newton, and G. L. Gross, "Mathematical 
Theory of Rocket Flight", McGraw Hill, New York (1947). 


% R. A, Rankin, Philosphical Transcations, Roy. Soc. of London, (A), 
241:457-585, (1949). 


"%4 P. Carriére, Memorial de l'Artillorio Française, 25:253-360 (1951), 
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automatic sensing and measuring control system, i.e., an optimalizing 
system which automatically holds the airplane at the measured optimum 
operating conditions. 

Of course, a skilled human operator controls the performance of 
a machine on the optimalizing principle: He watches the instrument 
readings of the inputs and outpute of the machine, and then uses his 
knowledge and experience to decide in what directions should the controls 
be adjusted. The adjusted inputs bring new output readingo which bave 
to be interpreted by the operator to determine whether the optimum 
operating condition is reached or exceeded. New adjustments of the 
control will have to be made. The continuous adjustment of inputs is 
the sensing process and the reading of the outputs is the feedback. How- 
ever, manually-controlled optimalizing systems are necessarily slow in 
response, and for complicated systems human skill, no matter how 
developed, is not sufficient. Automatic optimalizing control was con- 
ceived by C. S. Draper, Y. T. Li and H. Laning, Jr.* Its application 
to cruise control of airplane was discussed by J. R. Shull. 


min Principles of Optimalizing Control 


The heart of an optimalizing control system is the non-linear 
component which characterizes the optimum operating conditions. For 
simplicity of discussion, we shall assume that this basic component has 
a single input and a single output. For the time being, we shall neglect 
also any time effects and assume that the output is determined only by K .- nue 


vilu He input. Since there is an optimum operating point, output as a function 


of input has a maximum at & and KIM shown in Fig. 15.1. It is 
convenient to refer the output and the input to the optimum point and put 
the physical input as د‎ ag. and the physical output as ue . The 
optimum point is then the point x = 入 =0. The purpose of an optimal- 
izing control is then to search out this optimum point and to keep the 


* Y, T. Li, Instruments, 25:72-77, 190-193, 228, 324-327, 350-352 
(1952). C. $. Draper and Y. ۲۰ Li, "Principles of Optimalizing 


Control Systems and an Application to Internal Combustion Engine", 
ASME Publications (1951). 


+% J, R. Shull, Trans. I. R. E. (Electronic Computers), Dec. 1952, 
pp. 47-51. 
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m 
optimum filter design lay abanionimg the inadequate RC-circuits and 

hy actually using an analog computer or even a digital computer to 

serve as the filter. Then the theoretical optimum performance can 

be actually attained. However, tbe introduction of electro mechanical 
computer as a component of the filtering system certainly greatly 
increases the complexity of the overall system, and can be justified 

only in very critical cases. But if we have made the system very 
complicated at sime cost, we may ask whether we will actually obtain 

the very best performance. The optimum performance in the theory 
discusaed in the previous sections is only optimum within the limitations 
of the assumptions of the theory. For instance, two random signals 

with the same correlation function or the same power spectrum, 
according to the theory developed, require the same optimum filter. 
This is, in a sense, a certain looseness in design criteria, Surely, 

if we have more statistical information about the signal than just power 
spectrum, we should be able to distinguish these two signals and to 
improve our design by utilizing such additional knowledge. Then we 

can obtain even better performance than possible with the so called 
"optimum" filter. Itis evident that this generalized approach to the 
filtering problem must require more advanced theory of probability 

than we have used. The recently developed science of information مه‎ 
theory may also find important applications here. A beginning" has been 
made in this "probabilistic" approach to problem of detecting signal in 
noise. But much remains to be done. 

The Wiener-Kolmogoroff theory of optimum filter is based upon 
the mean square error criterion. By using this criterion, we essentially 
put the emphasis on minimizing the large errors without much consider- 
ation on the small errors. However in many occasions, we may be most 
interested in making the frequent error as small as possible, while not 
particular about making an infrequent large error. It is also possible 
that the probability function is very lopsided, with the mean far from 
the mode. For such cases, the mean square error criterion is entirely 


* See for instance, P. M. Woodward, I. L. Davies, Phil. Mag. 
41:1001-1017 (1950); Proc. I. R. E. 39:1521-1524 (1951); J.I, E. E. 
(London) 99(III); 37-51 (1952). T. G. Slattery, Proc. I. R. E. 
40:1232-1236 (1952). 


Chapter XVIII 


Control of Error 


In the preceding chapter, we have shown how the principle of 
ultrastability can make the control system insensitive to accidental 
errors and occasional failures of the components by the simple device 
of changing the characteristics of the system whenever instability occurs. 
Since an ultrastable system will automatically seek stability, the control 
system, when designed, actually embodies unstable fields of behavior as 
well às stable fields of behavior, In other words, during the design of 
an ultrastable system, we make no attempt to predict stability from 
instability, to separate the right fields of behavior from the wrong fields 
of behavior. Errors of behavior are merely treated as a probability, 
but otherwise unspecified. In this chapter, we shall approach the 
reliability of complex control system from a different point of view: We 
shall specifically introduce errors into the system and ask how should 
the system be designed so that the system will give satisfactory perfor- 
" mance in spite of the errors. That is, we wish to know how to control 
the error. 

U This subject of control of error is now in its early period of 
development. Only the control of error in most elementary operation 
can be discussed, And this is wholly due to J. von Neumann. * Our 
discussion in this chapter is then an exposition of Neumann's work. Its 
purpose is to serve as an introduction to this very important topic and 
to indicate the need for much further investigations. 


Reliability by Duplication. 

It is common knowledge that reliability of a system can generally 

3 be increased by the simple expedient of duplication, For instance, ifa 
simple system as shown in Fig. 18.1a has the characteristic that when 
it fails to operate, it merely gives no output. Then to guard against 

| probability of failure, we can duplicate the system with # identical 
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Eqs. (18. 27) and (18. 28) have a first term identical with Eq. (18. 1). 
The additional terms come from the imperfect elements and from the 
statistical distribution of errors. 

With any specified E, . 6 and m , Eqs. (18.26), (18.27), 
and (18.28) enable us to compute the distribution function of J , the 
fraction of activated output lines of the complete Scheffer stroke system. 
We can make this somewhat clearer by reverting to the notation of 
probability distribution functions. Thus for instance Eq. (18.26) is 
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The probability distribution function of 7 , W(Y;5,75%) , is thus the 


result of integrating with respect to g and A of the joint probability 
of E, /^, and لر‎ . Thus 
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We shall now show that under proper conditions we can obtain 
almost perfect performance of the multiplexed Scheffer stroke system 
by increasing 2 . Consider a given fiduciary level Š , Perfect por- 
formance means the implication of 75 Š , or non-activation of out- 
put, by F2/-d, 924-9 , or the activation of both inputs; the 
implication of JZ/,-S , by either کح"‎ „927-8 or F2/-d, 


